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STORAGE PROBLEMS WHEN DEMAND IS "ALL OR NOTHING" 



by 

D. P. Gaver and P. A. Jacobs 

Department of Operations Research 
Naval Postgraduate School 
Monterey, CA 93940 

0 . ABSTRACT 

An inventory of physical goods or storage space (in 
a communications system buffer, for instance) often experiences 
"all or nothing" demand; if a demand of random size D can be 
immediately and entirely filled from stock it is satisfied, 
but otherwise it vanishes. Probabilistic properties of the 
resulting inventory level are discussed analytically, both 
for the single buffer and for multiple buffer problems. 
Numerical results are presented. 



1 . INTRODUCTION 

The usual storage or inventory problems involve demands 
imagined to occur randomly, and to be capable of reducing any 
available stock to zero, or even beyond when backordering 
is permitted. Yet in many situations at least one component 
of total demand is "all or nothing; " that is, it reduces 
inventory only if it can be entirely satisfied by the inventory 
present, and otherwise seeks another supplier. Here are 
examples . 
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(a) A manufacturer's warehouse is filled with a certain item 

at the beginning of the selling season; let I denote 
the initial inventory. Suppose that demands occur as 
follows: a message is sent requesting that items 

be shipped from inventory, but only if the entire order 

can be filled. That is, the demand is satisfied if ^ / 

in which case inventory level is reduced to 1(1) = I - ; 

while if > I the inventory remains unchanged 

and 1(1) = I. Allowing for no replenishment, the second 
demand, of size D 2 t interacts with inventory 1(1), 
so that it is filled if D 2 £ 1(1)/ but is not placed if 
D 2 > 1(1). The process continues along these lines until 
the selling season is over and there are no more demands. 

(b) A buffer storage device used to contain messages prior to 
their batch transmission has capacity I. Messages of 
length {D^, i = 1,2,...} approach the buffer successively, 
and are admitted on an "all or nothing" basis, just as 

was true of demands for physical inventory in (a) above. 

Once again rejection will occur, and more frequently to 
large demands (messages) than to short ones. 

(c) A system of many buffer storage devices is used to contain 
messages prior to their batch transmission. Each buffer 
has capacity I. Messages of length {D^, i = 1,2,...} 
approach the device and are successively admitted to the 
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first buffer until there is a demand that exceeds its 



remaining capacity. The first buffer is left forever 
and the demand that exceeds the first buffer, plus 
successive demands, apply to the second buffer until one 
occurs that exceeds the remaining capacity. This demand 
then applies to the third buffer, and so on. As a 
result there will be some unused capacity in each 
buffer. 

In Section 2 we will discuss some models for the 
situations in examples (a) and (b) . We will compute such 
quantities as the distribution of the amount of inventory 
left at some time t and the distribution of the times of 
successive unsatisfied demands. 

In Section 3 we will consider a model for example c. 
We will derive equations for the limiting distribution of 
used capacity of a buffer and the expected used capacity of 
a buffer. It seems to be difficult to obtain simple analytic 
solutions to these equations, but we will present certain 
illustrative numerical results. 
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2. THE ONE-BUFFER INVENTORY PROBLEM 



Suppose that demands for available stock occur 
according to a compound Poisson process; if N^ is the ‘ 
number of demands that occur in (0,t], then {N^;t >^0} 
is a stationary Poisson process with rate A; the sizes of 
successive demands {D^} are independent with common 
distribution F. Assume that there are no replishments 
of inventory. Let t 0} denote the stochastic 

process describing available inventory at time t, and 
let {l (n) ; n = 0,1,...} be the stochastic process of 
available inventory following the nth demand. It is 
apparent from our assumptions that both {l(n)} 

are Markov processes. 

2.1. Functional Equations for the Amount of Available Inventory 
Let 

-si. 

(2.1) <|)(s,t)=E[e ] 



be the Laplace transform of the available inventory at time t. 
Similarly, let 



^^(s,n) = E[ . 



Properties of the available inventory can be studied in terms 
of <p and ij<; we begin by deriving an equation for <p. 
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observe that if one conditions on then in (t, t + dt) 

either no change in inventory occurs, an event of probability 
1 - XF(I^)dt + o (dt) , or a depletion of amoiint x occurs 
with probability A dt F(dx), x £ I^. Thus 

(2.2) E [exp (-sl^^^^) I 
-si 

= e [l-AF(I^) dt] + I exp [-S ~ x) ] Adt F (dx) + o (dt) 

0 

Now take expectations with respect to I^: 



(2.3) 4> (s , t + dt) 

= (j)(s,t)-A{E[e ^F(I^)]-E[/ exp{-s (I^-x)} F (dx) ]} dt + 



-si, 



After subtraction of (j>(s,t) from both sides, division by dt, 
and allowing dt 0 we find that <p must satisfy the 
equation 

(2.4) |1 = XE l^e ^ (e^^-1) F(dx) 

An analogous argument shows that 




(2.5) 



i|^(s, n+1) = \p(s,n) + E [e 



-si (n) 



/ 

0 



Kn) 



(e®^ - 1) F(dx) ]. 



o(dt). 
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Differentiation with respect to s at s = 0 , or a direct 
conditional probability argument, now produce equations for 
E [I^] and E [I (n) ] : 



(2.6) 


^ "t 

^ E[I ] = -XE[/ X F(dx) ] 

UX. t Q 


and 


Kn) 

E[I(n+l)] = E[I(n)] - E[/ X F(dx)] , 

0 



In general no explicit solutions for the expected 
values are available, but a simple upper bound results from 
rewriting (2.6) as follows. 



(2.7) 


^ "t 

jr E[I ] = -AE[I / ^ F(dx)] 

^ r 0 J-t 

> -AE[I^ ^(^t^ ^ 




> -AF(I) E[I^] , 



from which one sees that 

(2.8) E[I^] > Iexp[-AF(I) t] 

and similarly 

E [I (n) ] > I [1 - F(I) , 

so the expected available inventory declines by at most an 
exponential rate . 
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2.2. Explicit Solution When the Demand Distribution is Uniform 
Although Equation (2.4) seems to be quite intractable 
for most demand distributions , it can be solved completely 
when F is uniform: 

1 c > X 




and c > I. In this case (2.4) can be expressed as 




Ari-4>1.A 34) 

c L s J c 3s * 



In other words 4> satisfies a first-order (quasi) 
linear partial differential equation with initial condition 
<t)(s,0) = e . The usual procedure for solution, Sneddon 
(1957 ) r requires solution of two ordinary differential 
equations selected from among 



( 2 . 10 ) 



dt _ -ds _ d^ 

1 (A/c) “ ( A/c) [ (l-4>) /s] ' 



we find from the first and last two that 
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( 2 . 11 ) 



1 - 4 > _ 



s 




- c. 



so a general solution is given by 

(2.12) g^s + ^ t, ^ 3 = 0 » 



g being a function to be determined. Now the specified 
initial condition stipulated that at t = 0 



(2.13) 




/ 



SO 



(2.14) 




0 



which specifies (p at t = 0. But for t positive we 
replace s, the first argument of g at t = 0 by 
s + (X/c)t to obtain the solution 

(2 15) 1 - (j)(s,t) _ 1 - exp[-(s + (A/c)t)I] ^ q 

^ ' s ~ s + (A/c)t ' 

which gives the desired transform. Passage to the limit as 
s ->■ 0 in (2.15) shows that 
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(2.16) 



1 = 1 - exp [- (xt/c) I] 
(A/c) t 



This formula can also be derived by first finding an expression 
for the k th moment of and then employing a Taylor series 

argxament. 

In order to invert the transform in (2.15) note that 



(2.17) 



/ P{I^ 

0 ^ 



> x} dx 



1 - 4>(s,t) _ 1 - exp[-(s + (At/c))Il 
s s’ + ( a/c) t 



which is the transform of a truncated exponential distribution. 
Thus by the unicity theorem for Laplace transforms 



0 _< X < I, 

1 _< X . 

Note that the distribution of is absolutely continuous 

in the interval (0,1) but that there is a j\imp at I 
corresponding to the occurrence of no demand less than or 
equal to I in (0,tj: 



(2.18) P{I^ > x} = 




(2.19) 



P{I^ = 1} = expt- At(I/c)] . 
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2.3. The Expected Number of Satisfied Demands 

Supposing that an initial inventory, or storage capacity, 
I prevails, it is of interest to compute the probability that 
a demand is satisfied, and the expected number of demands 
satisfied in an interval of length t. First notice that if a 
demand of size D(t) appears at time t, at which moment 
is available, then 



is the conditional probability that the demand is satisfied. 
When F is uniform, as is presently true, we may remove the 
condition to find that 



If S(t) is the number of demands satisfied during the time 
interval (0,t], then since demands arrive according to a 
Poisson process with rate X, 



P{D(t) < = F(I^) 



P{D(t) ^ = E[F(I^)] = E 




1 - exp [- (At/c) I] 
At 



t 



t 



1 ~ exp [" (Au/c) I] 
Au 



(2.20) E[S(t)] = A / E[F(I^J]du = A / 




du 



= Y + + E (Ail) 

C 1 c 



where E^ ( • ) is an exponential integral; Abramowitz and 
Stegun (1965), and y = 0.5112... is Euler's constant. 
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2.4. The Time of the First Unsatisfied Demand and the Amount 



of Unused Inventory at That Time 

As before F is the common distribution function of 
the successive demands. Now let t be the time of the first 
unsatisfied demand. Then 

P{ T > t|N. = n} = P{D^ < I, D„ < I-Dw..., D <I-D, D J 

t 1— 2— 1' n— 1 n-1 

where F' ' denotes the nth convolution of F with itself. 

Hence 

(2.21) ■ P(t > t) = ^ e~^^ p(n) ^ 

n=0 

Explicit expressions for the distribution of t can 
be obtained in some cases. If F is uniform on [0,c] with 
c 2 then 

(2.22) P{T > t} = e"^^ Iq 




where Iq(z) is a modified Bessel function of the first kind 
of the zeroth order. In this case 



(2.23) 



E[t] = y exp{l/c} = y exp{l/2E[D]} . 
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If F is exponential with mean 1/y, then 



(2.24) 



P(t > t} 



I .1 e 



n 



-XI 



n=0 



n : 



k=n 



UliJl 

k; 



and 

(2.25) E[t] = i [1 + yl] = 1 [1 + . 

Note that if I is small relative to E[D], then the 
expected time to first unsatisfied demand when F is exponential 
will be greater than the expected time when F is uniform. How- 
ever for I large relative to E [D] the expected time for F 
exponential will be less than the expected time when F is 
uniform. 

Let be the amount of inventory present at the time 

of the n^^ unsatisfied demand. Then for 0 ^ i 

a 

P{Y > I - a} = / R(dy) F(I - y) 

0 

oo 

R(y) = I F^"^ (y) 
n=0 

F(I - y) = 1 - F(I - y) . 

Again explicit expressions for the distribution of Y^ 
can be obtained for some distributions F. If F is uniform 
on [0, c] for c > I, then 



(2.26) 

where 

(2.27) 
and 

(2.28) 
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(2.29) P{Yj^ ^ 1 - a] = 1 - - ) exp{^ a} . 



If F is a truncated exponential 



1 - e 



-yx 



(2.30) F(x) = 



1 - e 



-yl 



X < I , 



X > I , 



then 



(2.31) P{Yj^ >_ I-a} = 1 - [e"^^ - e”^^] [1 - e ^ exp{ ya [l-e“^^ ] } 



If F has an exponential distribution with mean 1/y , then 



(2.32) 



> I - a) = . 



In this last case the distribution function of Y can be 



n 



computed by induction quite easily and 



(2.33) 



P{Y > I - a) = . 



n — 



Hence when F is exponential 

(2.34) e[Y ] = ^ [1 - . 

n ny 

In principle similar results can be obtained for other distri- 
butions, but we have found no simple expressions. 
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2.5. Inventory Costs and Policies 

There are at least three monetary quantities which affect 
the profitability of an inventory policy over a fixed interval of 
time (0,t]: the selling price, p; the storage cost, a; and 

the cost of lost demands, b. If the storage cost a is 
charged just on the basis of I (something like warehouse size) 
then the total expected profit in (0,t] is 



for the case of uniformly distributed demands; see ((2.16) and 
(2.20)). One can nximerically find the maximum expected profit 
for this case; nothing explicit seems to be available. 



Z(I) = p(I - E[I^]) - al - bS(t) 




[1 - exp [- (At/c) I] ] 



- b{y + Jin (— I) + (^ I) } 

c 1 c 
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3. THE MANY-BUFFER STORAGE PROBLEM 



In this section we will study a model for the situation 
of example (c) in section 1. Messages are successively admitted 
to the n th buffer until there is a message length that exceeds the 
remaining capacity of the buffer. The total amount of this message 
is put in the (n+l)st buffer and the n th buffer is left forever. 
Successive messages are then put in the (n+l)st buffer until there 
is a message whose length exceeds the remaining capacity of the 
(n+1) st buffer; this message is put in the (n+2)nd buffer and 
so on. 

Let I denote the common capacity of the buffers and 
denote the length of message i. Assvune {D^} is a sequence 
of independent identically distributed random variables with 
distribution F having a density function f such that 
f(x) > d > 0 for X 6 [0,1]. Let R(x) = ^ (x) be 

the renewal function associated with F. If F(I) < 1, then 
we will assume that an incoming message to the currently used 
n th buffer of length greater than I is sent to the (n+l)st 
buffer; when it cannot fit into the (n+l)st buffer, then it is 
"banished," i.e. sent to some other set of buffers. The next 
message however will try to enter the (n+l)st buffer. If this 
message has length greater than I it is banished and the 
following message will try to enter the (n+1) st buffer; all 
messages of length exceeding I will be banished until one 
appears that is smaller than I and it will be the first 
entry in buffer (n+1) . 
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This model has been studied for demand distributions F 



with F(I) =1 by Coffman et al. (1978). Their approach was 

to study the Markov process describing the total amount of 

inventory or space consumed in successive buffers or bins. 

Here we study the process {L }, where L is the size of the 

n n 

demand that first exceeds the remaining capacity of the n th 
buffer; is a Markov process. Let 

K(x, [0,y]) = i ylL^ = x} . 



Note that 

P{L^ < y} = K(0, [0,y]) 

is the same as the sum of the forward and backward recurrence 
times at time I for a temporal renewal process with inter- 
renewal distribution F. Thus for y < I 



I 

(3.1) H^(y) = P{L^ 1 y> = / R(dz) [F(y) - F(I-z)]. 

I-y 

Note that for y < I 



(3.2) K(x, [0,y]) = - 



I-x 

/ R(dz) [F(y) -F(I-x-z) ] 

I-x-y 



if X < I-y; 



I-x 

/ R(dz) [F(y) -F(I-x-z) ] 

0 

if I-y < X < I; 



I 

/ R(dz) [F(y) -F(I-z) ] 

I-y 

if X > I . 
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Hence 



iR(I-x) - R(I-x-y)] F(dy) 



if X < I- 



IR(I-x) F(dy) 



if I > X 



(3.3) K(x,dy) = 



|R(y) F(dy) - / R(dz) f(y-z) + R(dy)F(y) 
0 

if X = I- 



[R(D - R(I-y) ] F(dy) 



if X > I 



Note that for some 0 < a < b < i, there exists a 6 > 0 
that for all x 

2 

K (x,dy) 2 <5 for y 6 [a,b] 

2 ^ 

where K (x,dy) = K(x,dz) K(z,dy). Hence hypothesis D 
on page 197 of Doob (1952) is satisfied. Thus, if 

K^(x,A) = ^ ^*^1 "" 

for all Borel subsets A, then 



(3.4) 



lim K^(x,A) = H(A) 
n -»■ » 



exists and further the convergence is geometric 

|K^(x/A) - H(A) I ^ 

for some positive constants a and y, y < 1 for all A 
Now let 



Y/ 

> i-y. 



y / 



such 
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H (x) = P{L e [0,x] |L = 0} . 
n n u 



Then a renewal argxainent can be used to show that for x I 



I 

(3.5) ^ / H^*R(dy) [F(x) - F(I-y)] 

I-x 

I 

+ [1-H (I)] / R(dy) [F(x) -F(I-y)] . 

" I-x 



Taking limits as n ~ it is seen that the distribution 
H(x) satisfies the following equation for x < I: 



I 

(3.6) H(x) = / H*R(dy) [F(x) - F(I-y)] 

I-x 

I 

+ [1-H(D] / R(dy) [F(x) - F(I-y)] . 
I-x 



Equations (3.1) and (3.6) can be simplified for certain 
specific distributions F. 



A. Exponential Demand s . 

For the exponential distribution with mean 1 and 
X £ I the equations are 



(3.7) Hj^(x) = 1 - e ^ - xe”^ 
and 

-X -X 

(3.8) H(x) = xe H(I) + H (x) - e / H(I-x+u)du . 

^ 0 
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B. Uniform Demands 



For the uniform distribution on [0,c] with c > I 
they simplify to 



(3.9) 

and 



H, (x) = expl^ (I-x) ] - (1 - ^) exp(^ I) 

J- ^ ^ c 



I-x 



(3.10) H(x) = i exp[^ (I-x)] / exp(- u) H(u)du 

c c Q c 



+ ^ H(I) - ^ exp(^ I) (1 - ^) / exp(- ^ u) H(u)du 



+ [1 - H(I) ] (x) , 



for X ^ I. Similar expressions hold for x > I, but they are 
unimportant in the present context. 

Equations (3.6), (3.8) and (3.10) do not seem to yield 

explicit answers. As a result we have solved (3.8) and (3.10) 
numerically by iteration using the system of equations 

_ X 

(3.11) f^n+1^^^ = xe ^ ® ^ / Hj^(I-x+u)du 

with as in (3.7) and 

1 1 I “X -1 

(3.12) ^^n+1 “ c / exp(- - u) (u) du 

c " c ®^P^c ■ c^ ^ exp(- I u) H^(u)du 

+ [1 - (I) ] (x) 

n 1 



with as in (3.9). 



For the cases carried out the convergence 
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is rapid; after n = 5 iterations very little change is noted 
and convergence has occurred, for most practical purposes. 

Next let be the amount of storage space used in 

the n th bin; the distribution of is denoted by 

G(x) = lim P{Y < x) = lim G (x) 

n — n 



is the long-run distribution. By probabilistic arguments 
and (3.4) 



(3.13) 



X 

G(x) = / H*R(dy) F(I-y) 
0 



X 

+ [1-H(D ] / R(dy) 

0 



F(I-y) 



where F(I-y) = 1 - F(I-y) and the long run average expected 
capacity of a bin that is actually used is 



I 

A = / X G(dx) . 
0 



For the case in which F is exponential with unit mean 



(3.14) A = I - [l-H(I)][l-e ^] - e ^ / e^ H(x)dx . 

0 



For the case in which F is uniform on [0,c] with c I 



I 1^1 

(3.15) A = -2 / H(u)d + exp(- I) / exp(- - u) H(u)du 

0 ^ 0 ^ 



+ H(I) [21 -cexp(- I) + c] + [-1 + c exp(^ I) - c] 
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Numerical solutions were obtained for equations (3.14) and 

(3.15) by first computing the probabilities H^(x), n = 1,2,..., 10 

iteratively from (3.7) and (3.11) for the exponential demand 

case, and from (3.9) and (3.12) for the case of uniform demands. 

Our technique was simply to discretize x: x^ = jh, h = I/N, 

N being the number of x-values at which H^(x) is evaluated 

(values of N from 200-1200 were utilized in order to obtain 

two-significant digit accuracy) . The integrals were then 

approximated by a summation, i.e. Simpson's rule. Having the 

values of H (x.) it is possible to calculate those of 
n 3 ^ 

H .1 (x.) , and from these the values of G (x) and the mean 
n+1 j n 

usage, E[Y^], may be calculated by numerical integration. In 
the case of exponential demand very simple upper and lower 
bounds were obtainable; such bounds were not tight enough to 
be useful for the uniform case. 

The following table summarizes the numerical results. 

We have compared demand distributions that result, as nearly 
as possible, in the same probability that an initial demand 
on an empty bin will be rejected. We have tabulated the 
expected level to which the bin is filled. It is interesting 

that the limited bin occupancy is 0.75 when a uniform demand 
over the range of the bin size is experienced. This result has 
been obtained analytically by Coffman et al. (1978); in that 
paper simple and elegant analytical expressions for G and H 
also appear for this case. The considerable similarity of 
the numbers in the rows of the table is notable; apparently the 
long-run bin occupancy is only slightly larger than is that 
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of the first bin, and the occupancy experienced for uniform 
demand is only slightly larger than for exponential. Further 
investigations to examine the reasons for this insensitivity 
would seem to be of interest. 

Expected Fraction of Bin Filled 

(f = E[Y ] ^ I) 
n n 

Rejection Probability Exponential Demand Uniform Demand 



F(I) 




fl 




f 

CO 


fl 


f 

oo 


0.00 




- 




- 


0.76 


0.75 


0.05 




0. 74 




0.75 


0.74 


0.74 


0.10 




0.69 




0.70 


0.72 


0.72 


0.15 




0.65 




0.66 


0.68 


0.69 


0.20 




0.60 




0.62 


0.64 


0.66 


0.25 




0.56 




0.58 


0.60 


0.62 
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